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$J 1 Abstract 

We derive the DGLAP and BFKL evolution equations in the N = 4 supersymmetric 
gauge theory in the next-to-leading approximation. The eigenvalue of the BFKL kernel 
in this model turns out to be an analytic function of the conformal spin \n\. Its analytic 
continuation to negative \n\ in the leading logarithmic approximation allows us to obtain 
residues of anomalous dimensions 7 of twist-2 operators in the non-physical points j = 
0, —1, ... from the BFKL equation in an agreement with their direct calculation from the 
DGLAP equation. Moreover, in the multi-color limit of the N = 4 model the BFKL and 
DGLAP dynamics in the leading logarithmic approximation is integrable for an arbitrary 
number of particles. In the next-to-leading approximation the holomorphic separability 
of the Pomeron hamiltonian is violated, but the corresponding Bethe-Salpeter kernel has 
the property of a hermitian separability. The main singularities of anomalous dimensions 
7 at j = —r obtained from the BFKL and DGLAP equations in the next-to-leading 
approximation coincide but our accuracy is not enough to verify an agreement for residues 
of subleading poles. 
PACS: 12.38. Bx 



1 Introduction 

The Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation [1, 2] and the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) equation [3] are used now for a theoretical description of 
structure functions of the deep-inelastic ep scattering (DIS) at small values of the Bjorken 
variable x. The structure functions are measured by the HI and ZEUS collaborations [4]. 
The higher order QCD corrections to the splitting kernels of the DGLAP equation are 
well known [5, 6]. On the other hand, the calculation of the next-to- leading order (NLO) 
corrections to the BFKL kernel was completed comparatively recently [7, 8, 9]. 

In supersymmetric gauge theories the structure of the BFKL and DGLAP equations 
is significantly simplified. In the case of an extended N = 4 SUSY the NLO corrections 
to the BFKL equation were calculated in ref. [9] for an arbitrary value of the conformal 
spin \n\ in a framework of the dimensional regularization (DREG) scheme. Below in the 
section 3 these results are presented in the dimensional reduction (DRED) scheme [10] 
which does not violate the supersymmetry. The analyticity of the eigenvalue of the BFKL 
kernel as a function of the conformal spin \n\ gives a possibility to relate the DGLAP and 
BFKL equations in this model, as we show below (see also [9]). Further, in both schemes 
(DREG and DRED) this eigenvalue has the property of the hermitian separability which 
is similar to the holomorphic separability (see Section 3). 

Let us introduce the unintegrated parton distributions (UnPD) ip a (x,k 2 _) (hereafter 
a = q,g,<p for the spinor, vector and scalar particles, respectively) and the (integrated) 
parton distributions (PD) f a (x ) Q 2 ) 



In DIS Q 2 = —q 2 and x = Q 2 /(2pq) are the Bjorken variables, k± is the transverse 
component of the parton momentum and q and p are the virtual photon and hadron 
momenta, respectively. 

The DGLAP equation has a clear probabilistic interpretation and relates the parton 
densities with different values of Q 2 : 



where the second term in the r.h.s. of the equation is the Mellin convolution of the 
transition probabilities Wb^ a (x) and PD fb(x, Q 2 ). Usually the first and second terms are 
unified with the use of the splitting kernel Wb^ a [3]: 



It is well known, that equation (2) is simplified after its Mellin transformation to the 
Lorentz spin j representation: 





(1) 
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fa(j,Q 2 )=J2^(j)f b (j,Q 2 ), 



dlnQ 2 
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where 



fa(j,Q 2 )= f dx X J- l f a (x,Q 2 

Jo 



are the Mellin momenta of parton distributions. The Mellin moment of the splitting 
kernel 



lab(j) = I d,X X 3 l W b -> a {x) 

Jo 



coincides with the anomalous dimension matrix for the twist- 2 operators 1 . These oper- 
ators are constructed as bilinear combinations of the fields which describe corresponding 
partons (see Eq. (3) below). 

On the other hand, the BFKL equation relates the unintegrated gluon distributions 
with small values of the Bjorken variable x: 

d r 

ip g (x,k 2 ± ) = 2u(-k 2 ± ) ip g (x, k]_) + / d 2 k' ± K(k±, k' ± )(p g (x, k 2 ± ), 



d\n(l/x) 

where uj{— kj_) < is the gluon Regge trajectory [1]. 

Let us introduce the local gauge-invariant twist-two operators: 
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(3) 

where the spinor \1> and field tensor G Pili describe gluinos and gluons, respectively. The last 
expression is constructed from the covariant derivatives D p of the scalar field $ appearing 
in extended supersymmetric models. The symbol S implies a symmetrization of the tensor 
in the Lorentz indices /ii, fj,j and a subtraction of its traces. 

The matrix elements of Of, ,, and O" „ are related to the moments of the distri- 
butions f a (x, Q 2 ) and Af a (x, Q 2 ) for unpolarized and polarized partons in a hadron h in 
the following way 

f^dxx^Ux.Q 2 ) = <h\n^...n^O; i ^.\h>, a=(q,g,<p), 
f Q dxx 3 ' l AUx,Q 2 ) = <h\n^...n^O; u ^.\h>, a = (q,g). (4) 

Here the vector n M is light-like: n 2 = 0. In the deep-inelastic ep scattering we have 



1 As in Rcf. [7, 9], the anomalous dimensions differ from those used usually in the description of DIS 
by a factor (-1/2), i.e. lab (J) = (-l/2) 7 f 6 ' s (j). 
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The conformal spin \n\ and the quantity 1 + uo (u is an eigenvalue of the BFKL ker- 
nel) coincide respectively with total numbers of the transversal and longitudinal Lorentz 
indices of the tensor 0" „ with the rank 

J =l+u + \n\. (5) 
Namely, we can introduce the following projectors of this tensor 



where the complex transverse vector l± is given below 

11 = + ll = l ± p = l ± q = 0. 

It is important, that the anomalous dimension matrices jabij) an d %b{j) for the twist- 
2 operators O a ^ l 4l . and 0^ l 4l . do not depend on various projections of indices due to 
the Lorentz invariance. But generally for the mixed projections the Lorentz spin j of the 
tensor is less than the total number of its Lorentz indices J. 

The matrix elements of the light-cone projections are expressed through solutions f a 
of the DGLAP equation (see Eqs.(4)). On the other hand the mixed projections 

n' 11 ...^^...!^ < p \°%, ..., N \ p > ~ £ dxx" J d 2 k ± (J^J <p g (x,k 2 ± ), (7) 

can be found at small u> in terms of solutions ip g of the BFKL equation. 

Thus, it looks reasonable to extract some additional information concerning the par- 
ton ^-distributions satisfying the DGLAP equation from the analogous /^-distributions 
satisfying the BFKL equation. Moreover, due to the fact, that in an extended N = 4 
SUSY the /^-function vanishes, the 4-dimensional conformal invariance may allow us to 
relate the Regge and Bjorken asymptotics of scattering amplitudes. 

Our paper is organized as follows 2 . In Section 2 we discuss the relation between 
DGLAP and BFKL equations in the N = 4 model obtained by an analytic continuation 
of the eigenvalue tv°(n, v) in \n\. In Sections 3 and 4 we review shortly our previous 
results for the BFKL kernel [9], rewrite them in a framework of the DRED scheme and 
investigate their properties. Section 5 is devoted to an independent calculation of anoma- 
lous dimensions in the leading and next-to-leading approximations for N=4 SUSY with 
the use of the renormalization group. In Section 6 we discuss the relation between the 
DGLAP and BFKL equations in the NLO approximation for this model. Appendices A, 
B and C are devoted to the derivation of some results used below. A summary is given 
in Conclusion. 

2 Some results were presented on the XXXV Winter School (see [11]) 
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2 Anomalous dimensions of twist-2 operators and their 
singularities 

In the leading logarithmic approximation fermions and scalars do not give any contribution 
to the BFKL equation and therefore its integral kernel is the same at all supersymmetric 
gauge theories. Due to the Mobius invariance in the impact parameter space Jf the 
solution of the homogeneous BFKL equation has the form (see [12]) 

(\m/ „. \m 
_PV2_\ 
P10P20/ VP10P20/ 

where 

m = l/2 + iu + n/2, m = 1/2 + iv - n/2 

are the conformal weights simply related to eigenvalues of the Casimir operators of the 
Mobius group. We use in -Ey, n (pio, pio) the complex variables pk = %k + Wk in the trans- 
verse subspace and the notation p ki — p k — Pi- 

For a principal series of the unitary representations the quantities v and n are respec- 
tively real and integer numbers. The projection n of the conformal spin \n\ can be positive 
or negative, but the eigenvalue of the BFKL equation in LLA [2] 

u = u\n,v) = ^(*(i)-ite<lf + + M)) (8) 

depends only on |n|. The Mobius invariance takes place also for the Schrodinger equation 
describing composite states of several reggeized gluons [13]. 

It is important, that the above expression has the property of the holomorphic sepa- 
rability [14] 



u/'(//>) +u°(m) , uj°(m) = ^-^(2^(1) - *(m) - *(1 - m)) (9) 



,0/ 



8tt 2 



due to the relation 



*(m) + ^(1 -m) = *(m) + - m) . 

Moreover, the more general property H = h + h* , [h,h*]—0 is valid also for the Hamil- 
tonian of an arbitrary number of reggeized gluons in the multi-colour QCD N c — > 00 [14] . 
For the case of the Odderon constructed from three gluons the holomorphic hamiltonian 
h is integrable [15]. Furthermore, in the multi-color QCD the BFKL dynamics for an 
arbitrary number of reggeized gluons is also completely integrable [16]. It is related to the 
fact, that in this case the holomorphic Hamiltonian h coincides with the local Hamiltonian 
for an integrable Heisenberg spin model [17]. The theory turns out to be invariant under 
the duality transformation [18]. Presumably the remarkable mathematical properties of 
the reggeon dynamics in LLA are consequences of the extended N = 4 supersymmetry 
[19], which is one of the reasons for our investigation of the BFKL and DGLAP equations 
in this model. It was argued in Ref. [19] that generalized DGLAP equations for matrix 
elements of quasi-partonic operators [20] for iV = 4 SUSY are also integrable at large 
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N c . For the case of the Odderon a non-trivial integral of motion discovered in Ref. [15] 
was used to construct the energy spectrum and wave functions of the three-gluon state 
[21]. A new Odderon solution with the intercept j equal to 1 was found in [22]. Recently 
the energy spectrum for reggeon composite states in the multi-color QCD was obtained 
through the solution of the Baxter equation (see refs. [23] and [24]). The effective action 
for reggeon interactions is known [25]. It gives a possibility to calculate next-to-leading 
corrections for the gluon trajectory and reggeon couplings [7]-[9]. 

The solution of the inhomogeneous BFKL equation in the LLA approximation can be 
written as the four-point Green function of a two-dimensional field theory 

. x , , ,_ K,\n\(p~w,P2^i)E* M (pro,p^o) 



/OO r 
dvC{v, \n\) / d 2 p - 
-oo J 



uj — u°{\n\, v) 



where C(v, \n\) is expressed through an inhomogeneous term of the equation with the use 
of a completeness condition for E Vi \ n \ (see [12]). 
For pit — > p2> we obtain por ~ ~Pv2> and therefore 

< m) m) 40$) > ~ E r d "C(v, \n\) EvM{ ^'^ pVvpfv (10) 

n J -co to — uj y\n\^v) 

jp I >■ >\ / \ M/2 

£>v w ,\n\{PlV , P2V ) | ,l+2iv u / PV2' \ 



n 



c \y^ M) — T7^\ — \ — \pi'2> 



u'(\n\,v u ) \P\>2>) 
where is a solution of the algebraic equation 

uj = u; (|n|, v) 

with Im Vu < 0. 

The above asymptotics has a simple interpretation in terms of the Wilson operator- 
product expansion 

where 

1 

r w = - + iv u 

is the transverse dimension of the operator Vu ^ n \{p~y) calculated in units of the squared 
mass. This operator is a mixed projection 

of a gauge-invariant tensor O^,...^ with J = 1 + u + |n| indices. Note, that because 
T w is real in the deep-inelastic regime pu — > the operator Vui! \ n \(p~i') belongs to an 
exceptional series of unitary representations of the Mobius group (see [26]). 

The anomalous dimension 'y(j) obtained from the BFKL equation in LLA (see (8)) 
has the poles 
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r w = i + !^-7(j), ^')Uo = f^- (ii) 

Note, that there are also solutions of the BFKL equation corresponding to operators 
with the shifted anomalous dimensions 

r w (A;) = V u + k (A; = 1,2,3,...). 

They contain the covariant laplacians D 2 and have higher twists. 

The canonical contribution 1 + \n\/2 to the transverse dimension corresponds to 
the local operator 

in accordance with the fact, that in the light-cone gauge A^n^ = the tensor 

has the transverse dimension equal to 1. 

The local operator 0„ w ,| n | for |n| = 1, 2,... has the twist higher than 2 because its 
anomalous dimension 7 is singular at uj — * 0. Indeed, for integer \n\ > such singularity 
is impossible for twist-2 operators, because in this case for uj — » the total number 
of their indices 1 + uj + |n| would coincide with the physical value of the Lorentz spin 
j — 1 + |n| > 2. Moreover, the twist of Vuj \ n \ should grow at large |n|. For the case 
uj > \n\ > this operator is a mixed projection of the tensor 

C C C D DC 

In particular, for n = its twist equals 2. For uj < \n\ the corresponding operator is 

C C 1 C D D D C 

" p£t 1 " CTi • • • ^ ^ cr^ _ 1 ^ (7 W ^ a w + 1 • • • 1J a„ ^ pp 1 + w 

and its anomalous dimension has a singularity at uj — > 0. 

In the previous paper [9] we suggested to use an analytic continuation of the BFKL 
anomalous dimension 7(|n|, uj) to the points \n\ = — r — 1 (r = 0, 1, 2, ...) to calculate the 
singularities of the anomalous dimension of the twist-2 operators at negative integer j 

j = 1 +uj + \ n \ -> -r . (12) 

Because for fixed |n| 7^ — r — 1 and u; — > the quantity 7(]n|,u;) describes higher twist 
operators, to obtain the anomalous dimension 7 for the twist-2 operators one should push 
\n\ to — r — 1 sufficiently rapidly at uj — > 

A(|n|) = \n\ + r + 1 = Ci(r) cj 2 + <5 2 (r) uj 3 + ... . (13) 

Note, that for iV = 4 SUSY there are several multiplicatively renormalizable twist-2 
operators (3) constructed from bosonic and fermionic fields. Their anomalous dimensions 
can be obtained from the universal function 7 Mm (j) (see Eq. (71)) by a shift of its 
argument j — > j + k with k — 0, 1, 2, 3, 4. It is related to the fact, that these operators 
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enter in the same supermultiplet. We assume in accordance with Ref. [9], that from the 
BFKL equation one is able to calculate the singularities of 7 nm (j). 

In LLA one can easily derive from the BFKL eigenvalue tu°(n, v) in the limit (12) the 
following result for ^(j) 

l(j)\^-r = (14) 

for all r = —1,0,1,.... Further, for the BFKL equation in LLA the fermions are not 
important, but generally they give non- vanishing contributions to the residues of the 
poles for 7(7) in the DGLAP equation even in LLA. Therefore the result (14) for the 
anomalous dimension 7 um (j) can be derived only for a certain generalization of QCD. 

In our previous paper [9] it was shown, that only in an extended N=4 supersymmetric 
Yang-Mills theory the anomalous dimension, calculated in the next-to-leading approxi- 
mation from the BFKL equation, can be analytically continued to negative \n\. Therefore 
the above result for 7(7) in LLA should be valid only for N = 4 SUSY. Indeed, using the 
conservation of the energy-momentum stress tensor (corresponding to the condition 
7(2) = 0) to fix a subtraction constant in the expansion of 7(7) over the poles at j = —r, 
we obtain : 

2 A T 

70') = f^7 LLA (j), 1 LLA U) = 4(¥(1) - ¥(j - 1)) 

in an agreement with the direct calculation of 7 um (j) in this theory (see [19, 28] and Sec- 
tion 5 below). It is important that this anomalous dimension is the same for all twist- 2 
operators entering in the N = 4 supermultiplet up to a shift of its argument by an inte- 
ger number, because this property leads to an integrability of the evolution equation for 
matrix elements of quasi-partonic operators [20] in the multi-colour limit N c — > 00 (see 
[19]). 

In the NLO approximation there is a more complicated situation. Namely, the anoma- 
lous dimension 7 um (j) has the multiple poles A7 ~ a 2 (j + r)~ 3 at even r, which is related 
to an appearance of the double-logarithmic corrections ~ (a In 2 s) n in the Regge limit 
s — > 00 at upper orders n of the perturbation theory. 

The origin of the double-logarithmic terms can be understood in a simple way using 
as an example the process of the forward annihilation of the e + e~ pair in the pair 
in QED [27]. Let us write the scattering amplitude for this process in the following form 

_i_ j_ 2 
/e+e— = 47ra em 7f7 ' la A(s, k 2 ± ), a em = ^ , s = -2p e k , 

S 47T 

where p e and —k are momenta of the incoming electron and positron respectively and 
k± is the transverse part of k. A(s, k±) is the Bethe-Salpeter amplitude for the electron 
interaction with the positron having the large virtuality k\ ^> m 2 . The integral equation 
for this amplitude in the double-logarithmic approximation (a ern \n 2 s ~ 1 , a em <^ 1) can 
be presented as follows [27] 

A(s, kl) = 1 + ^ £ ^ \l ^A(s>, k'l) , (s = min{ S , sk'l/kl}) , (15) 
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where we took into account, that the double- logarithmic contribution appears from the 
integration region 

k' 2 k 2 
s' < s , — p > — , s' > k'l > m 2 . 
s' s 

One can search its solution in a form of the double Mellin transformation 
Putting this anzatz in the equation we obtain a simple expression for f w (7) 



3 



= (,-2 7 ),^ 

, , <-*em J- 1 J- 



Z7T 



7 7 



The pole of fu{l) in 7 is situated at 



•^H 1 -^'- (17) 

Expanding this expression for the anomalous dimension 7 in the series over a em we shall 
generate the triple pole A7 ~ al m /uj 3 . 

It is important to note, that the pole contributions 

1 1 
+ - 

uj — 7 7 

in the denominator of ^(7) are similar to singularities of the expression 

w~2tf(l)-tf(l + |n|+w-7)-tf(7) (18) 

(at \n\ = —1, —3, ...) appearing for uj — > in the eigenvalue of the BFKL kernel of N = 4 
SUSY in a modified Born approximation (see Eqs. (44) and (45) below). For odd negative 
j the double-logarithmic contributions are absent because the kernel of the Bethe-Salpeter 
equation has an odd number of factors k± leading to a cancellation of one logarithm after 
the integration of / (Pk^/k 2 ^ over the azimuthal angle tp. 

It is obvious from the above QED example, that the Bethe-Salpeter equation within 
the double-logarithmic accuracy can be derived both from the BFKL-like and DGLAP- 
like equations by imposing more accurate constraints on the region of integration. In 
the first case one should restrict k' 2 by s' from above, which results in the substitution 
—7 — > uj — 7 in the argument of the first \l/-function (see the above expression for uj ). In 
the second case one should restrict s' by s from above, which results in the substitution 
7 _1 — > 7 _1 + (uj — 7) -1 and leads to a presence of the triple poles A7 ~ a 2 /uj 3 in the 
two-loop anomalous dimension. This procedure is similar to the rapidity veto approach 
invented by Bo Andersson and other authors to explain the appearance of the large NLO 
corrections to the BFKL kernel [29, 30]. 



3 The expression for / w (7) in the published version of this paper (see Nucl.Phys. B661 (2003) 19) 
was not right because it did not satisfy all boundary conditions contained in the integral equation (15). 
The correct expression (16) is derived in the Appendix D of this hep-ph version and is sent by us as 
an Erratum to Nucl.Phys. B. We thank J. Bartels and M. Lyublinsky who drew our attention to this 
problem. 
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3 NLO corrections to the BFKL kernel in N = 4 
SUSY 

To begin with, we review shortly the results of Ref. [9], where NLO corrections to the 
BFKL integral kernel at t — were calculated in the case of QCD and supersymmetric 
gauge theories. Only N = 4 SUSY is discussed in details. 



3.1 The set of eigenvalues 



The eigenvalues of the homogeneous BFKL equation for the N = 4 supersymmetric gauge 
theory can be written as follows [9] 



uj = 4a 
where 7 = \ — iv and 



X(n, 7) + -Xin, 7) + S(n, 7) J a 



X(n, 7 ) = 2tt(l)-tt(M)-*(l-M), 



(19) 



(20) 



with 



0(71,7) = ^ , (M) + M/"(l-M) + 6C(3)-2C(2) X (n,7) 
- 2$(|n|, 7 )-2$(|n|,l- 7 ) 



\ n \ N 

M = 7 + — , M = 7- — 

[ ^ 2 ' f 2 



(21) 



Here ^(z), *&'(z) and ^"(z) are the Euler \1/ -function and its derivatives, respectively; 
a = g 2 N c /(lQir 2 ) is expressed in terms of the coupling constant g in the DREG scheme. 
The function $(|n|,7) is given below 4 



„M-1 



$(H, 7 ) = -/ o -(*'(B^)-C(2))+Li 2 (-x)+Li 2 (x) 

+ ln(x)U(|n| + l)-*(l) + ln(l + x) + Vf^l 
V fc=1 fc + \n\ 



+ E 



E 

fc=0 



k ^(k + \n\) 2 
-l) k+1 



-If) 



k + M 



V'(k + \n\ + 1) - ty'(k + 1) + (-l) k+1 [p'(k + \n\ + 1) + (5'{k + 1)) 



(22) 



4 Notice that the representation of 3>(|n|,7) in terms of the sum over k in ref. [9] contained a misprint: 
the factors (— l) k+1 in the last sum of (22) were substituted by (— l) k . 
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and 



Note, that the term 



1 

4 



-E 



:-v, 



r =o + r ) 



(23) 



gX(",7) 



appears as a result of the use of the DREG scheme (see eq. (29) in [9]). It is well known, 
that DREG scheme violates SUSY. The DREG procedure could be used provided that 
the numbers of boson and fermion degrees of freedom would coincide (see [31]). But in 
its standard version the number of degrees of freedom for vector fields is 4 + 2e 5 and does 
not depend on e for other particles. Usually (see discussions in [32]- [34]) it is possible to 
restore the equality between boson and fermion degrees of freedom in SUSY by adding 2i 
additional scalar fields. Formally it is equivalent to the dimensional reduction (for N = 4 
SUSY) from the space-time dimension 10 to the dimension A + 2i and we call it below the 
DRED scheme. The change of the number of scalars from 6 to 6 — 2e does not violate the 
cancellation of the singular contribution ~ e^ 1 to the coupling constant (the /^-function is 
zero for N = 4 SUSY in the DRED scheme 6 ) but leads to an additional contribution to 
uj of the order of a 2 . Indeed, the correct result in the framework of the DRED scheme can 
be obtained from eq.(19) by a finite renormalization of the coupling constant (see also G. 
Altarelli et al. and G.A. Shuler et al. in [34]) 



_ J-—2 

a = a + -a , 
o 



(24) 



which eliminates the term ~ a 2 x in (19). For the new coupling constant a in the DRED 
scheme the above expression for u can be written in the following simple form 



4a 



X(n,l) + <J(n,7)a 



3.2 Hermitian separability of the Bethe-Salpeter kernel 

Following the method of Ref. [28] and using the results of the previous section we can write 
the above NLO correction to the BFKL equation in the form having the property similar 
to the holomorphic separability (cf. [14]). This property can be called the hermitian 
separability, because it guarantees the reality of the eigenvalue u for real v. 

Indeed, using the results of Appendix A we can obtain for the most complicated 
contribution to 5(n, 7) 



$(|n|, 7) + H\n\, 1 - 7) = X(n, 7) (?{M) + (3'(1 - M 
+$ 2 (M) - p'(M) [^(1) - *(M)] + $2(1 — M) — p'(l - M) - m(l - M) 



5 Here £ = (D — 4)/2 and D is the space-time dimension. 

6 In the DRED scheme the /3-function vanishing was demonstrated at first three orders of perturbation 
theory (see [33, 34] and references therein). 
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where x(. n il) is given by Eq.(20). 

Thus, we can rewrite the NLO corrections S(n, 7) in a (generalized) hermitian separable 
form (providing that Uo is substituted by u, which is valid with our accuracy): 

5(n, 7 ) = 0(M) + 0(l-M)-g(p(M)+p(l-M)), (25) 

ujq = 4a(2#(l)-tf(M)-tf(l-M)) (26) 

and 

p(M) = (3\M) + ic(2) , (27) 

0(M) = 3C(3) + V(M) -2$ 2 (M) +2/3'(M)(*(l) - (M)) . (28) 

In a simpler way the hermitian separability takes place in the eigenvalue equation for the 
corresponding Bethe-Salpeter kernel 

1 = ^(2*(l)-*(M)-*(l-M)+a(0(M)+0(l-M))) 

-2a (p(M) + p(l - M)). (29) 

In the right hand side the first contribution corresponds to its singularity at I — — 1 in 
the Gribov-Froissart representation generated by a pole of the Legendre function Qi(z) at 
uj = 1 + I — > and the last term appears from the regular part of the Born contribution. 
Note, that M and 1 — M coincide with the anomalous dimensions appearing in the 
asymptotic expressions for the BFKL kernel in the limits when the gluon virtualities are 
large: q\ — > 00 and q\ — > 00, respectively. Because 1—M = M*, the hermitian separability 
guarantees the symmetry of a; for the principal series of the unitary representations of the 
Mobius group (with m — 1/2 + iv + n/2) to the substitution v — > —v and the hermicity 
of the BFKL hamiltonian. 



3.3 The violation of a generalized holomorphical separability 

Eqs.(25)-(28) show the hermitian separability with two contributions depending on M = 
7+ \n\/2 and M* = 1 — M — 1 — 7+ \n\/2. Although the above expressions are symmetric 
to M <-> M*, it is not clear, if we can symmetrize 5(n, 7) to the substitution \n\ <-> — |n| 
as it was done in the Born approximation to obtain a true holomorphic separability 
(see Eq.(9)). Using properties of polygamma-functions and Eqs. (A4) and (A5) one 
can attempt to present Eqs.(25)-(28) in a separable form symmetric to the substitution 
m <-> m (i.e., respectively, M <-> M for both positive and negative n). From the results 
of Appendix B the NLO correction 8(n, 7) can be written in the form 

5(n, 7) = Mm) + 0(m) — — ° ( p(m) + p(m) 

2a V 

o /cos 2 (m7r) cos 2 (m7r)\ r(2 )/ x /0 ^ x 
\sm a (m7r) sur(m7r) / 
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where 



4%, 7) = + <-!)") -»<*?) 

sm [Mn) 1 



2 1 ^(m) cos(m7r) fy(m) cos(m7r) 



(1 + (-1)") e(n) 7r v 7^ : 7 - r , 2 J y (31) 
\ sm 2 (m7r) snr(m7r) / 



with 



e(ra) = 



+1, n>0 
-1, n<0 



The term 5^ (n, 7) violates the holomorphic separability. Because the contribution 
proportional to ujq in (30) depends not only on M and 1 — M we have the situation simi- 
lar to the case of quantum anomalies. Namely, the property of the hermitian separability 

of 8(n, 7) is responsible for the violation of the holomorphic separability in 5^\n, 7). The 
(2) 

anomalous term 62 \ n il) is zero f° r °dd n, where S(n, 7) coincides with its analytic con- 
tinuation to corresponding negative \n\. Note, that for the Pomeron only even values of 
n are physical due to the Bose symmetry of its two-gluon wave function. The colorless 
composite state of three reggeized gluons with the /-coupling in the color space and odd 
n has an anti-symmetric wave function and can give a large contribution to the small- a; 
behaviour of the structure function g2{x) (see [36]). 



3.4 Asymptotics of "cross-sections" at s — > 00 

As an application of obtained results we consider the cross-section for the inclusive pro- 
duction of two pairs of particles with their mass ji in the polarized 77 collisions (see 
[2, 9]): ' ' ( ^ 

= "emO 2 -^— (a (s) + (cos 2 ??- -)<r 2 (s)), 

where a em is the electromagnetic fine structure constant, the coefficient cr (s) is propor- 
tional to the cross-section for the interaction of unpolarized photons and a 2 (s) describes 
the spin correlation depending on the azimuthal angle d between the polarization vectors 
of colliding photons in their cm system. 

The asymptotic behavior of the cross-sections (T^{s) (k = 0, 2) at s — > 00 corresponds 
to an unmoving singularity of the t-channel partial wave f w {t) ~ (u — uJk) 1 ^ 2 situated at 



Uk = 4a 



X [k, - + a5 [k, - 



9tt 5 / 2 



1/2 



<T 2 (S) 



32 V / 7C(3) (\n(s/s )) 
9.32 v /7C(3)-8(ln( s / So )) 1/2 



= 4ax \k, - 



l + 0(a 



1 — ac ( k, - 

^ 2 



l + 0(a 



(32) 
(33) 
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Here the symbol 0(a) denotes unknown next-to- leading corrections to the impact factors. 
Note, that strictly speaking for the N = 4 supersymmetric gauge theory we can not 
consider electro-magnetic interactions without the SUSY violation and it could be more 
natural to investigate the interaction of gravitons or their superpartners. But here we 
want only to illustrate the relative magnitudes of the radiative corrections to the BFKL 
equation in QCD and in N=4 SUSY. 

Using our results (25) and (26), we obtain in the N = 4 case the following numerical 
values for X {k, 1/2) and c(k, 1/2) (k = 0, 2): 



X 



4 



4 In 2, 



X 



4(m2-l), 



(34) 



2C(2) + 



= 2C(2) + 



21n2 



HC(3)-32Ls 3 (|)-147rC(2) 



7.5812, 



(35) 



2(ln2 - 1) 



11C(3) + 32Ls 3 (^) + 14ttC(2) -32 In 2 



6.0348 , 



(36) 



where (see [35, 37]) 



Ls 3 (z) = -^ln 2 2sin(|) 



dy. 



Note, that the function Ls 3 (a;) appears also in contributions of some massive diagrams 
(see, for example, the recent papers [38] and references therein). 

The LLA results (34) coincide with ones obtained in ref. [2]. As it was shown in 
[7, 9], in the framework of QCD the NLO correction c QCD (0, 1/2) is large and leads to 
a strong reduction of the value of the Pomeron intercept (see recent analyses [39, 29, 
41, 30] of various resummations of the large NLO terms). Contrary to c® CD (0, 1/2), the 
correction c(0, 1/2) is not large (c QCD (0, l/2)/c(0, 1/2) « 3.5). Because in N=4 SUSY 
the /9-function is zero, it is natural to interpret the large correction to the intercept of the 
BFKL Pomeron in QCD as an effect related to the coupling constant renormalization. 
It seems to support the results of refs. [39, 40] (see also the recent review [30] and 
references therein) concerning a large value for the argument of the running QCD coupling 
constant at high energy processes. The corrections c® CD (2, 1/2) (see [7, 9]) and c(2, 1/2) 
for N = 4 are small and do not change significantly the LO value [2] of the angle-dependent 
contribution. 



4 Non- symmetric choice of the energy normalization 



Analogously to refs. [7, 9] one can calculate eigenvalues of the BFKL kernel in the 
case of a non-symmetric choice for the energy normalization so in eq.(15) related to the 
interpretation of the NLO corrections in the framework of the renormalization group (cf. 
[7, 9]). For the scale so = q 2 natural for the deep- inelastic scattering we have respectively 
in DREG-scheme 



LO 



4a 



X{n, 7 - tJ) + Q*( n > 7) + S(n, 7)) a 
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= 4a 



X(n, 7) + ( 7) + <J(n, 7) J a 



in DRED-scheme 



w = 4a 



X(n,7) + 5(n,7) ,a 



where 



and 



5(71,7) = <J(n,7) -2x(n,7)x'(n>7) ( 37 ) 



x '(n, 7) = ^X(n, 7) = -*'(M) + - M) . 

Note, that here 7 does not coincide with the anomalous dimension of the higher-twist 
operators with |n| > 1 (see below). 

4.1 Limit 7^0 for n = 

By considering the limit 7 — > of the BFKL eigenvalue one can obtain for n = (see also 
the analysis in [9]) 

X(0,7) = - + 0( 7 2 ), 
7 

1 ~ -qDREG 

-X(n, 7 ) + 5(0,7) = + C + 0( 7 2 ), 



7 

tjDRED 

5(0,7) = + C + 0( 7 2 ), 

7 



(38) 



where 



B 



DREG 



= 3' 5 



DRED 



= 0, C = 2C(3). 



(39) 



According to refs.[7, 9] with the use of eqs.(38) and (39 ) one can calculate the anoma- 
lous dimensions 7 of the twist-2 operators at u — > (i.e. near j = 1) respectively 
in DREG-scheme 



7 = 4 a 
in DRED-scheme 

7 = 4 a 



(i + 0( U )) +B ( 



j^DREG 



+ 0(1) + a 



(i + 0(w) ) + a (^ + 0(1) ) + a »(£ + („-.)) 



C 



(40) 



(41) 



Thus, in the framework of the DRED scheme the singular contribution ~ a 2 jto to the 
anomalous dimension is zero. These results will be used below for the calculation of NLO 
corrections to 7 from the DGLAP equation (see Eqs. (68) and (71)). 
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4.2 Hermitian separability for a non- symmetric normalization 
and the symmetry 7 <-> J — 7 

For the scale So = q 2 the expression for u as a function of the anomalous dimension 
7 = 1/2 + iv + \n\/2 (correctly defined at general n (see (11)) can be written in the 
following form 



^ = 4a (x fw,7- -y) + ^ ( n ' 7 ~ 2))' 



where x( n ?7) an d S(n, 7) (for the non-symmetric choice of so) are given by Eqs. (20), 
(21) and (37), respectively. 

By summing double-logarithmic terms in all orders of the perturbation theory (see 
discussion in Section 2) we can write uj in the "Lorentz invariant" form 



UJ 



4a(W(l) - #(7) - + \n\ + u - 7) + A( n, 7) a) , (42) 



where the Lorentz spin j of the corresponding operators (for negative \n\) is given by (12). 
For n = the Lorentz spin of the twist- 2 operator is j = 1 + uj. 
In the above expression 



A(n, 7 ) = «jU,7-i^J +2 



* , ( 7 ) + * / (|n| + I- 7) 



xU-Y-y ■ (43) 



Using the analysis of the subsection 3.2 one can present uj as follows 

uj = 4a (2 (1) - (7) - *(|n| + 1 + - 7) + e) , (44) 
where e is written for uj — > in the form 

e = + M -7)) +0(0(7) + 0(1 + N -7)), 

= 4a(2*(l) - #(7) -*(|n| + l -7)) + 0(a 2 ) . (45) 

Here 

P(7) = *'(7)-p(7) = *'(7) - ^(7) - = 2E (7+ 1 2fe)2 -^C(2) 

and 0(7) is given by Eq.(28). 

In accordance with the hermitian separability of the BFKL kernel established in the 
subsection 3.2 the eigenvalue equation for the corresponding Bethe-Salpeter equation can 
be written as follows 

An 

1 = _(2*(l)-*( 7 )-*(J- 7 ) + a(0( 7 ) + 0(J-7))) 
u 

+2a(p( 7 ) + p (J- 7 )) , J = 1 + \ n \ +u, (46) 

where J is the total number of tensor indices of the local operator which does not coin- 
cide generally with its Lorentz spin j. Note that without the shift l + |n|— 7 ^J — 7 
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of the argument of \I/-function at the LLA level, the symmetry of the eigenvalue to the 
substitution 7 — > J — 7 is violated by the term x'( n ? 7 ~ M/2)- 

Let us calculate A(n, 7) near its singularities at small 7. To begin with, we consider 
7 — > for the physical integer values |n| > 0: 

A(n, 7) -> + +^(c(n) + 2*'(|n| + l)), (47) 

where 

c (n) = *'(|n| + 1) - *'(1) - + 1) + /3'(1) . (48) 

Therefore by solving the equation u = iv(n, 7) one can obtain 

7 = f (l + " ( *(1) - *(M + 1))) + (*(!) " *(W + !) + | C W) • ( 49 ) 

At n = the correction ~ a to 7 is absent, but for other n we have the large correction 

A7 = 4a (^(1) - + 1)) , 

having the poles at 1 + uj + \n\ — > — r, which leads to a contribution changing even the 
singularities of the Born term. The explanation of this effect is related to the presence 
of the double-logarithmic terms A7 ~ ci 2 /uj 3 near the points j = 0,2, ... (see Sections 2 
and 6). We remind, that for positive integer |n| we calculate the anomalous dimensions 
7 of the higher twist operators (with an anti-symmetrization between n transversal and 
1 + iv longitudinal indices). The singularities of the anomalous dimensions of the twist- 2 
operators can be obtained only in the limit when A(|n|) tends to zero more rapidly than 
(v (see Eq. (13). 

5 Anomalous dimension matrix in the N = 4 SUSY 

The DGLAP evolution equation for the moments of parton distributions for N = 4 SUSY 
has the form 

^2lab(j)fb(j,Q 2 ) (a,b = q,g,<p), (50) 
k 

J2% b (j)&MQ 2 ) (a,b = q,g), (51) 
k 

where the anomalous dimension matrices 7 a b(j) and 7 a b(j) can be obtained in the form 
of expansions over the coupling constant a 

% b (j) = a.^(j) + a*.^(j) + 0(a 3 ). (52) 



a 

dlnQ 2 
d 



fa(j,Q 2 ) 



d\nQ< 



A/ a (j,Q 2 ) 
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In the following subsections we will present the results of exact calculations for the 
leading order coefficients 7^(j) and 7^(7) an d construct the anomalous dimensions of 
the multiplicatively renormalizable twist-2 operators. In the NLO approximation the 
corresponding coefficients 7„^(j) and 7^(j) were unknown (see, however, [42]). It is 
important, that the form of the LLA anomalous dimension matrix of the multiplicatively 
renormalizable operators in N = 4 SUSY is rather simple because the result is expressed in 
terms of one function. Taking into account this universality related to the superconformal 
invariance, the existing information about the anomalous dimensions of twist-2 operators 
in the QCD case, the known NLO corrections to the BFKL kernel and an experience 
in integrating certain types of the Feynman diagrams (see, for example, [43, 44]), we 
derive below the expressions for the NLO anomalous dimensions of the multiplicatively 
renormalizable operators in N = 4 SUSY. This result is checked by direct calculations of 
the matrix elements 7^(7) an d t|^(j) [42]- 



5.1 LLA results for the anomalous dimension matrix in N = 4 
SUSY 

The elements of the LLA anomalous dimension matrix in the N = A SUSY have the 
following form (see [28]): 



for tensor twist-2 operators 



2 1 1 



3 j + lj' >mKJ ' j + 1' 

1 I 

J-l 3. 



7$(j) = 4 (*(1) - *(j + 1)) , 7 $(j) = 4 ( — - 7 ) , (53) 



for the pseudo-tensor operators: 



2 



J + 1 3 J 

rj 0) (3) = sM + AV = 2 /2 



= 4{*(l)-*0- + l) + -^---l. (54) 



3 3 + 1/ V./ J + 1, 

1 L 



J + 1 J 



Note, that in the = 4 SUSY multiplet there are also twist-2 operators with fermion 
quantum numbers but their anomalous dimensions coincide up to an integer shift of the 
argument with the above expressions for the bosonic components (cf. ref. [20]). 
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5.2 Anomalous dimensions and twist-2 operators with a multi- 
plicative renormalization 

It is possible to construct five independent twist-two operators with a multiplicative renor- 
malization. The corresponding parton distributions and their LLA anomalous dimensions 
have the form (see [28]): 

flti) = fgU) + fqU) + UU) ^ fq,g,V^ 

7f(j) = 4(*(l)-*(j-l)) = -45i(j-2)= 7 f (J), (55) 

fn(j) = -2V-i)f a U) + f q V) + lu + i)fM~f! a ,M 

ifhj) = 4(*(l)-*(j + l)) = -45i(j)=7J 0) (j), (56) 



fill (j) = - fg U )+fqU) — U U ) ~ fq,g,<p U ) > 



-'—Im+m- 3 -^- 1 

ifkj) = 4(*(l)-*( J + 3)) = -4S 1 (j + 2)= 7 i%'), (57) 



f IV (j) = 2Af B V) + Af q V)~Af+ g V), 

ifiij) = 4(*(l)-*(j')) = -45i(j-l) = 7f0'), (58) 

Mi) = -(j-l)Af g (j) + J -±^Af q (j)^Af- g (j), 

l { v\j) = 4(*(l)-*(i + 2))^-45 1 (i + l)^7i° ) 0'), (59) 

Thus, we have one supermultiplet of operators with the same anomalous dimension 
l LLA {j) proportional to ^(1) — ^(j — 1) 7 ■ The momenta of the corresponding linear 
combinations of parton distributions can be obtained from the above expressions fk(j) 
by an appropriate shift of their argument j in accordance with the corresponding shift of 
the argument of 7fc(j). Moreover, the coefficients in these linear combinations for N = 4 
SUSY can be found from the super-conformal invariance (cf. Ref [20]) and should be the 
same for all orders of the perturbation theory in an appropriate renormalization scheme. 

The momenta of three multiplicatively renormalizable twist-2 operators for the unpo- 
larized case are 

/jv(j') = a g f g (j) + a q f q (j) + a v f v (j) 

where the coefficients can be extracted from above expressions (55)- (59). If we insert 
this anzatz in the DGLAP equations (50) the following representations for the correspond- 
ing anomalous dimensions 

7^0) = iffo) + %ffo) + %Su) 

Ug Ug 



7 A similar result for another type of operators was obtained in ref. [45] . 



18 



(lq (lq 

can be obtained. Eqs.(60) lead to relations among the anomalous dimension matrix 
lab U) ( a ? b = g,q, tp) which should be valid also in the NLO approximation up to effects 
of breaking the superconformal invariance (see [42] and references therein). 

Analogously for the set of two multiplicatively renormalizable operators in the polar- 
ized case 

Af N (j) = a g Af g (j)+a q Af q (j) 
we can derive the following relations 

fflV) = lf 9 \j) + ^lf 9 \j) = lf q \j) + H° q \j) ■ (61) 

(Ig (lq 

So, we have nine equations for the matrix elements in the case of the usual partonic 
distributions and four equations for the polarized distributions, which determines com- 
pletely the anomalous dimension matrices 7^0) (a, b = g,q,ip) and 7^(i) {a, b = g,q) 
in terms of their eigenvalues in LLA 

7 i 0) (j) = -45x(j =F 2), 7S 0) (J) = -45x0-), l ( ±\j) = -45x(j =F 1) • 
This procedure is considered in details in Appendix C. 



5.3 NLO anomalous dimensions and twist-two operators with 
the multiplicative renormalization 

We have the following initial information to predict the NLO anomalous dimensions of 
twist-two operators with the multiplicative renormalization in N=4 SUSY. 

1. As it was shown in the previous subsections, the LLA anomalous dimensions are 
meromorphic functions having the poles at j = — r, r = —1,0,1,.... Moreover, there is 
only one basic anomalous dimension r ) LLA (j) and all others can be obtained as ^ LLA (j+m), 
where m is an integer number. It is useful to choose (see Eq. (14)): 

1 LLA U) = 4(tt(l) - *V - 1)) = -4Si(j - 2). 

Then, r y LLA (j) has a pole at j — > 1 and vanishes at j = 2. One should keep the above 
universality and linear relations among the matrix elements also for the NLO anoma- 
lous dimensions ^^{j) and 7„^(j) because it is a consequence of the super-conformal 
invariance. It means, that we should construct only the basic NLO anomalous dimension 

i NLO U)- 
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2. There are known results for the NLO corrections to the QCD anomalous dimen- 
sions. 



3. In the MS'-scheme with the coupling constant a and in the MS'-like scheme with 
the coupling constant a (i.e. in the scheme based on DRED procedure), the terms ~ ((2) 
should disappear in the final result for the forward Compton scattering (see [46] -[48]). 
Therefore the terms ~ £(2) are cancelled at even j in anomalous dimensions for the struc- 
ture functions F 2 and F L (related to the unpolarized parton distributions) and at odd j in 
anomalous dimensions for structure functions g\ and F 3 (related to the polarized parton 
distributions). 

4. From the BFKL equation in the framework of DRED scheme (see (20), (21) and 
(22)) we know, that there is no mixing among the functions of different transcendentality 
levels % 8 , i.e. all special functions at the NLO correction contain the sums of the terms 
~ 1/n 1 (i = 3). More precisely, if one will introduce the transcendentality level of the 
functions in accordance with the complexity of the terms in the corresponding sums 

# ~ 1/n, ~ 0' ~ C(2) ~ 1/n 2 , ~ (3" ~ C(3) ~ 1/n 3 , 

then for the BFKL equation in LLA and in NLO the corresponding levels are % — 1 and 
i — 3, respectively. 

Because in N=4 SUSY there is a relation between the BFKL and DGLAP equations, 
the similar properties are assumed to be valid for anomalous dimensions themselves, i.e. 
the basic functions r ) LLA (j) and "i NLO (j) should be of the types ~ \jf with the levels 
% — 1 and % — 3, respectively. The only exception could be for the terms appearing in the 
Born approximation, because such contribution can be removed by an approximate finite 
renormalization of the coupling constant. The LLA basic anomalous dimension is given 
above in terms of S\{j — 2). Then, the NLO basic anomalous dimension 7 Ari °(j) can be 
expressed through the functions: 

S±i(j-2), S ±k , t (j-2), ((k)S ±l (j-2), C(<) 

(here i — 3 and k + I = i), where 

S 2 ,iU) = ti%)- (62) 

m=l 171 

S-2U) = E - (i - (-i) J ')k(2), 

m=l Z 

8 Note that similar arguments were used also in [44] to obtain analytic results for contributions of some 
complicated massive Feynman diagrams. The method was based on a direct calculation of several terms 
in the series over an inverse mass with taking into account its basic structure found earlier in [43, 44] by 
considering few special diagrams with the use of the differential equation method [49]. 



20 



j I i \m o 

S-sU) = £ ^ - (i - (-1) V (3) ' 

m=l 1,1 * 

S-2,iU) = H) J EM-5iH-(1-(-1)^C(3) (63) 

rn=l 171 

The functions S^U), S-s(j) and S -2,1(3), which were introduced recently in [50], 
coincide (up to an opposite sign) with the functions K 2 (j), K 3 (j) and if 2,1(7) considered 
in [47, 51] and used in our previous paper [11], i.e. S-i(j) = — iQ(?), S- 2 ,i(j) = — if 2,1(7')- 
Note, however, that our definition (63) of S-i(j) and S_k,i(j) coincides with the functions 
of [50] only at even values of j. Expressions (63) after the shift of the summation index 
m — > m — j can be analytically continuated from even to complex values of j, which 
reproduces j3'(j + 1) (23) and associated functions. 

Note that the terms ~ £(2) should be absent for even values of j and for odd values 
of j in the unpolarized and in polarized cases, respectively, in accordance with item 3. 

Moreover, there is an important observation: the function 

S-iU) = (-1)'' £^/--(l-(-l)')ln2 

m=l 171 

does not contribute to the QCD anomalous dimensions and the Wilson coefficient func- 
tions (see [5] and [47, 51], respectively). Further, as it was shown in [47], the terms 
~ S-i(j) cancel in the final results for the diagrams describing the longitudinal Wilson 
coefficient function. 

5. The terms 

S ±l (j-2)/(j±m) k (j + l = i) (64) 

are absent in the anomalous dimension 7 Ari °(j). There are two reasons for this conclusion. 

Firstly, these terms have additional poles in the points j = =Fm. But such poles 
should cancel, if we start with the BFKL equation and obtain , y NLO (j) by the analytic 
continuation to |n| = — 1 — r. Indeed, using such procedure one can not obtain the 
doublets of poles. 

The second reason comes from the consideration of the multiplicatively renormalizable 
linear combinations (58) and (59). If, for example, in the polarized case functions (64) 
contribute, then we shall have the terms ~ (j ± m) _1 in one combination and the terms 
~ (j ±m±2) _1 in another combination. However, from the direct calculation of the NLO 
anomalous dimensions in the polarized case (see [6]) we know that only the terms ~ j -1 
and ~ (j + l) -1 are present in these combinations. 

So, terms (64) should be absent in the universal NLO anomalous dimension in the 
N=4 SUSY case. 

Note, that the sums S'_ 2 (j), S , _ 3 (j) and S , _ 2 .i(j) (see Eqs.(63)) appear explicitly in 
calculations only at even values of j in the unpolarized case (and/or at odd values of j, 
after the replacement (— iy — > (— in the polarized case). The analytic continuation 
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to complex values of j for the functions S±i(j) (i = 1, 2, 3) and S-2.i(j) in Eqs. (62) and 
(63) can be done easily (see [47, 52]) by a replacement of the sums Z)m=i i n the r - n - s 
with the difference J2m=i ~ Em=i+j' They are expressed in terms of the polygamma- and 
associated functions: 



SiU) = + !)-*(!), 
W = 



-2,1 



(J) 



(i-1)! 



+ 1) - 



(— l) i-1 1T ,i_l 



(i-1)! 



* ,_1 (j + l) + C(0 (<>1), 



2C0) - £ ^ = 2C0) - E ' 



00' 
(-1) 



m =i (m + j) 2 ~ ~! (m + j) 2 

l)-/3(l)=/5(j + l) + ln2, 



+ !)-*(!) 



(65) 



i-l r 



(i-1)! 



0.-1^ + 1)_ ^-1(1 



(-1 



^(j + ij + CH) (<>i), 



(i-1)! 



g (-ir>Si(m + j + l) _ 5_ C(3) 



m=0 



E 



(m + j + l) 2 

^(m + j + 2) - *(1) 



m=0 (™ + i + !) 2 
where C(-r) = (2 x - r - l)C(r) (r > 1), C(-l)=h2. 



^C(3) = 3T'0- + 1)-|C(3),(66) 



6. Further, the NLO anomalous dimension 7 7VLO (j) is equal to a combination 
of the most complicated contributions (i.e. the functions with a maximal value of the 
transcendentality level % — 3) for the QCD anomalous dimensions (with the SUSY relation 
for the QCD color factors Cp — Ca = N c ). 

These most complicated contributions (with % — 3) are the same for all QCD anoma- 
lous dimensions (coinciding in N=4 SUSY) [5, 6] (only the NLO scalar-scalar anomalous 
dimension is not known yet): 



l$ QCD {j) 



7. 



(l)QCD 



gg 



U) 



1 1 



qq 



U) 



16( g ) A;' QU 2) • 



1 { £ QCD U) 
{C F = C a = N c ) , 



where we omit less complicated contributions and 



QU) 



2 ( SsU) + S-sU) ) + ( 5 2 (j) + 5- 2 (j) ) - 5- 2> i(j) • 



(67) 



(68) 



Using Eqs. (65) and (66), the function Q(j) can be rewritten in terms of the polygamma- 
and associated functions: 

QU) = lc(3)-P"U + i) + \(*"U + i) + P''U + i) 

+ (nj + 1) - m) (^c(2) - *'u + 1) - + 1)) 



(69) 



= lc(3)-P"U + i) + ^" 

+ i(*0' + l)-*(l))(c(2)-*'( 



J + l 
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Thus, for N = 4 SUSY the NLO universal anomalous dimension , y NLO {j) has the form 

7 ™(i) = 16Q(.7-2) (70) 



7. We could add the term ~ ((3) to the r.h.s. of (68), but due to the condition 
7 ArLO (j = 2) = it cancels. 

So, the universal anomalous dimension 7(7) in two first orders of the perturbation 
theory for N=4 SUSY is 



where 7 iLj4 (j) and 7 Ari °(j) are given by Eqs. (14) and (70), respectively. All other 
anomalous dimensions can be obtained as 7 Mm (j + m), where m is an integer number. 

Thus, the above arguments allow us to construct the NLO corrections to anomalous 
dimensions in N=4 SUSY, which were unknown earlier. We check these results by direct 
calculations [42] and reproduce the anomalous dimension 7 ArLO (j). Note, however, that in 
[42] the coupling constant a (see Eq.(24)) is used, which is responsable for an appearence of 
the additional contribution l/3 7 LO (j) to the above obtained NLO anomalous dimension 
(70). This contribution can be absorbed in a by a finite renormalization of the coupling 
constant (see the subsection 3.1). Note also that in the dimensional reqularization scheme 
the eigenvalues of this matrix are not expressed only in terms od the function 7 um (j) with 
the shift j — > j + m of its argument. 



Using our knowledge of the anomalous dimensions we can construct the solution of the 
DGLAP equation in the Mellin moment space in the framework of N=4 SUSY. 

A. Polarized case 

The polarized parton distributions are splitted in two contributions: 



l um U) 



ai LLA (j) + a^ NLO U), 



(71) 



5.4 DGLAP evolution 



A/,, 9 (j,Q 2 ) = A/+(j,Q 2 ) + A/- 9 (j,g 2 ), 



(72) 



where 



at LO 




:(o) 




(73) 



at NLO 





(74) 
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where 



/ ~(i) - \ 

~(1) . 

+ MZ 9 NLO {j,Ql). (75) 

Here the anomalous dimensions 7±^ and 7^. are related to 7^ (a, b — q,g) as follows: 



5?>0) 7«0) J I «0) t!?(;) 1 1 (76) 



a a -1 

with V and U given in eq.(C4). 

Notice that only anomalous dimensions 7±± are important for N = 4 SUSY at the 
order 0(a 2 ) because they contribute to the (^-evolution of parton distributions. 

The anomalous dimensions 7^, which were not calculated in the previous section, can 
be found from eq.(76) provided that the gluino-gluon polarized anomalous dimensions 
labii) ( a ib = AS are knowm The anomalous dimensions 7±^ give contributions at 
the order 0(a 2 ) only to the normalization factors A/^ WLO (j, Qg)- They appear in the 
Q 2 -dependent part of Af^ g (j,Q 2 ) at the level 0(a 3 ) in the following form: 



~(0) ~(0) 

7± -It 



B. Unpolarized case 

The unpolarized parton distributions are splitted in the three parts: 



where 
at LO 



U(i,Q 2 ) = E fl 9 Jj>Q% (77) 

»=+ -,o 



f q , 9 Jj,Q 2 ) = f^QDl^Y (7f= 45!(j=F2), 7? ) = 4S 1 (j)), (78) 

at NLO 



(0) 
2 X 7- 'a 



fL^ Q 2 ) = fiZFV, Ql) ' ' (7E = 16Q(j =F 2), 7$ = 16Q(j)) (79) 

As in the previous case A, only the anomalous dimensions 7^ are important at the 
level 0(a 2 ) in N=4 SUSY. The anomalous dimensions 7^ (i 7^ k) contribute at 0(a 2 ) 



2x T 
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level only to the normalization factors (j, Qo)- They are related to the qluino-gluon- 

scalar anomalous dimensions 7^0) {a, b — g,q, tp) as follows: 



\ 7?!(J) 



7o-0) 



7^) \ 



l { 99 \j) 



I g<] \J ) I qq^ \J ) I <pq 
V 7gip \3) Iqip \3) ^y<p<p 



7g(i) 

7 W (J) 



where V and V -1 are given by eq.(ClO). The anomalous dimensions 7^ (i 7^ A;) can 
be obtained from the above equation provided that the gluino-gluon-scalar anomalous 
dimension matrix 7^(j) (a, b — g,q, tp) is known. 



6 Relation between the DGLAP and BFKL equa- 
tions 

As we pointed out in subsection 3.6, in the case of iV = 4 SUSY the BFKL results (20) and 
(21) are analytic in \n\ and one can continue the eigenvalues to the negative values of \n\. It 
gives a possibility to find the singular contributions to anomalous dimensions of the twist-2 
operators not only at j = 1 but also at other integer non-physical points j = 0, —1, —2.... 
In the Born approximation for the anomalous dimension of the supermultiplet of the 
twist-2 operators we obtain r y unt = 4a (^(1) — ^(j — 1)) which coincides with the result 
of the direct calculations (see [19, 28] and the discussions in subsection 4.2 ). Thus, in 
the case of iV = 4 the BFKL equation presumably contains the information sufficient for 
restoring the kernel of the DGLAP equation. Below we investigate the relation between 
these equations in the NLO approximation. 



6.1 DGLAP approach 

Let us start with an investigation of singularities of the anomalous dimensions of twist-2 
operators which were obtained in a direct way at the previous Section. 

By presenting the Lorentz spin j as uo — r, where r = —1, 0, 1, ... and pushing uo — > 
one can calculate the singular behavior of the universal anomalous dimension 7 um (j) (71). 
Note, that in our discussion of the BFKL equation we used the more general definition 
j = 1 + uo + \n\ for the rank of the Lorentz tensor O^..^. (see (12)), where 1 + 00 is 
the number of its longitudinal indices and \n\ is the conformal spin. Therefore strictly 
speaking in all expressions discussed in this subsection uj should be substituted by j + r. 

For the special functions contributed to 7 LLA (j) and r y NLO (j) at j = uj — r — > — r, r > 
we have the following expansion in uj: 

(_l)i C (/ + l)_ 5l+1 ( r + l) 

(-l)'C(J + 2) + S I+2 (r + l) 

Ss(j-2) = -^ + S 3 (r + 1) 



S 1 (j-2) = — + S 1 (r + 1)-J2^ 1 



S 2 (j-2) = -_-<? 2 (r + l)-£(Z + l)u;< 
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2^ o u 



1=1 



(-l)'C(Z + 3) - Si+3(r + 1) 



(80) 



S-2U ~ 2) 



S-s(j ~ 2) 



= (-!) r+1 {^ + C(-2) (1 - (-ID + 5_ 2 (r + 1) 

OO r -| n 

+ EC + iV (-l)'C(-(i + 2))+S-( W )(r+l) 
= ("l) r+1 R + C(-3) (1 - (-1)0 - S_ 3 (r + 1) 



+ E 



(1 + !)(* + 2), , 



UJ 



1=1 



(-l)'C(-(Z + 3))-5_ (I+3) (r + l)]} 



(J " 2) 



(-ir +i (- 



la; . 
5 



C(2)-5_ 2 (r + l) 

<(3) + §C(3) (1 - (-l) r ) - S- 2>1 (r + 1) " 5_ s (r + 1)} . (81) 
Presenting our results for r ) LLA (j) and 7 7VLO (j) in terms of these functions, we obtain 



LLA/ 



7~ 0) = 4 
7 ™O0 = 8 



- - S x (r + 1) - S 2 (r + 1)uj + 0(uj 2 

UJ 

+ (-1)1 _ 2Sl(r + !)(! + (-!)') 



(82) 



o;- 



u;^ 



((1 + (-1DC(2) + 2(-l)'S 2 (r + 1)) - + O(o;°) 
N ' u; 
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X i -2S 1 {r + l)±-S 2 {r + l)i + 0{^ 
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S 2 (r + lH + 0(u; ) 



if r = 2A; 

if r = 2Jfe + 1 , 



(83) 



where 5* 2 (r) = £(2) + S 2 (r). So, the double-logarithmic poles ~ uj~ 3 appear in the case of 
even values of r (see discussions in Section 2). 

Note that the functions S-i(r + 1) and S- 2 (r + 1) do not contribute to 7 L (j) in the 
limit j — > — r, r > 0. The absence of S-i(r + 1) is explained by the fact, that the quantity 
S-i(j — 2) does not appear in the basic NLO anomalous dimension (see discussion in the 
subsection 5.3, item 4). 

On the other hand, the absence of the term S_ 2 (r + 1) is related to the following 
important property of 7. As it follows from eq.(81), the functions S- 2 (j — 2), S-s(j — 2) 
and S-2,i{j — 2) giving contributions to 7 Ari °(j) (see eqs. (70) and (68)) are responsible 
for the different asymptotics of it at j — > — r, r > for even and odd r (see the r.h.s. of 
eq.(83)). But the combination of these functions S- 2 (j — 2), S- 3 (j — 2) and S- 2 ,i(j — 2) 
contributes to 7 7VLO (j) in such way, that the function S_ 2 (r + 1) is absent in the r.h.s. of 
eq.(83). 

Thus, for 7 Mm (j) we obtain at j = uj — r — > — r (r > 0) 9 



l un \3) = G 



- - S x {r + 1) - u S 2 (r + 1) + O{uo 2 ) 

UJ 



9 Because really the expansion parameter is 4a (see, for example, Eqs. (82) and (83), we shall consider 
below G — 4a as a coupling constant. 
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+ G 



2 J £ - 25x(r + 1)£ - S 2 (r + 1)1 + O( W ) if r = 2k 



S 2 (r + 1)^ + 0^) 



if r = 2k + 1 



(84) 



Note that 7 nm (j) at j = u — r — > — r can be presented as a solution of the following 
Bethe-Salpeter equations 



1 = G 
for r = 2/c > and 

1 = G 



1 - uj S^r + 1) - to 2 S 2 (r + 1) 
7(w - 7) 

1 ^(r+lj+w^r + l) 



+ 1) 



(85) 



7 



+ S 2 (r + 1) 



+ ... 



(86) 



for r = 2/c + 1 > 0. 

As it was shown in the subsection 5.4 the anomalous dimensions, responsible for the 
Q 2 evolution of parton distributions in the framework of N = 4 SUSY model, have the 
following form 

in unpolarized case 



7±(j) 
7o (j) 



G7i 0) 0') + G 2 7i 1 i(i) + ... = 1 um (j + 2 T 2) 
G^°\j) + G 2 ^(j) + ... = 7 um (j + 2) 



and in polarized case 

l±{j) = G^\j) + G 2 ^l(j) + ... = 7™(J + 1T2). 



J7) 



(88) 



Thus, for unpolarized and polarized cases, the corresponding anomalous dimensions 
have the double-logarithmic poles u~ 3 at even and odd r, respectively. The arguments of 
regular terms and of the functions in front of poles uj~ % are shifted by an integer number. 
By chosing certain values of this number we obtain various 7 s (j). 



6.2 BFKL approach 

Let us investigate a possibility to obtain the residues of the anomalous dimension 7 urw (j) 
in the poles at j = —r (r = 0, 1,2, ... ) from the BFKL equation. Its eigenvalue u can 
be analytically continued to the negative points \n\ = — r — 1. One can present u as a 
solution of the equation (see (44)) 

to = G (2 (1) - (7) - *(|n| + 1 + u - 7) + e) , 

having the property of the hermitian separability corresponding to the symmetry between 
7 and |n| + 1 + to — 7 (see (46)). The expression for e is given in Eq. (45). 

To begin with, we note that providing that j = \n\ + 1 + uo (see Eq. (12)) the 
r.h.s. of the BFKL equation contains the ^-function with its argument equal to j at 
7^0, whereas r j LLA {j) obtained from the DGLAP equation contains \l/(j — 1). The 
functions contributing to the NLO correction have an analogous shift of their argument 
in comparison with the functions appearing in the NLO corrections to 7 um (j) (84). 
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In the Born approximation ^ LLA ~ G the reason for the difference of arguments of 
^-functions appearing in the DGLAP and BFKL approaches can be easily understood. 
Starting from the BFKL equation one reproduces the universal anomalous dimension 
7™ of twist-two operators by an analytic continuation of its eigenvalue uj to the points 
j — \n\ + 1 + to — to — r (r = 0, 1, ...). The residue of the pole of r y um at j — 1 is 
calculated also in accordance with the DGLAP equation but without any continuation 
(see (40), (41)). In the BFKL approach the term regular at j — > — r does not contain 
the contribution from the pole at j = 1 and equals —GSi(r) instead of —GSi(r + 1). It 
is important, that the anomalous dimension obtained from the BFKL equation depends 
on two parameters uj and |n| and therefore we can impose a certain relation (see (13)) 
between two asymptotics uj — > and \n\ + r + 1 — > 0. Another ambiguity is related to the 
existence of double-logarithmic terms near j = —r (see below). As a result, in particular 
in LLA we can obtain an agreement between the regular contribution to r y unt and its pole 
singularities. 

In the next-to-leading approximation the situation is similar but not so transparent. 
As in the case of LLA to obtain the residues of the poles for 7 um (j) in the BFKL approach, 
we continue the eigenvalue uj to j = — r (r = 0,1,...). But this continuation does not 
contain the regular contributions from the singularity at j — 1. Another problem is that 
to calculate the residues of the non-leading poles G 2 uj~ 2 and G 2 uj~ 1 at uj ~ j + r — > for 
7 um one should continue the BFKL equation from the region of its applicapability w < 7 
to the double-logarithmic region 7 ~ uj. The property of the hermitian separability (46) 
of the BFKL kernel is helpful for this purpose, because it relates the small parameters in 
one combination uj + \n\ — r — 1 — 7. It turns out, that in the intermediate region 7~wwe 
should rearrange the BFKL equation in such way that the Born contribution will contain 
the double-logarithmic terms. This resummation procedure leads to a renormalization 
of the next-to-leading correction ~ G 2 . Namely, to avoid the double-counting we should 
subtruct from this correction the sum of poles containing the square of the most singu- 
lar Born contribution. It turns out, that the subtracted terms are small in the BFKL 
region w < 7, but they are essential for the non-leading poles G 2 uj~ 2 and G 2 uj~ 1 of the 
anomalous dimension. These subtructions can be found by calculating the next-to-leading 
contribution to the BFKL kernel with a greater accuracy in comparison with ref. [7], but 
we leave this problem for future publications. Below we discuss in more details the origin 
of above ambiguities. 

Let us write down the eigenvalue of the BFKL kernel near the singularity at 7 = for 
the physical conformal spins |n| = 0,2, ... (see Eqs. (42), (43) and (47)) 

- = G (i - »(1 + M) + »(!)) + G 2 ( *(D-*(1 + M) + fM) . (fa) 

/,From this equation one can obtain easily the anomalous dimensions (49) at integer 
positive |n|. The analytic continuation of uj to the negative integer points \n\ — > — r — 1 
contains divergencies due to the following formulae 

(1) - *(|n| + 1) -> i,, 1 - - 5i(r) - (M + r + 1) S 2 (r) , 
\n + r + 1 
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c(n) - V~ ( ^ 2 + S 2 (r) + (-l)^_ 2 (r) - 1 (1 + (-1)') C(2) , 

where ^(r) was introduced after Eq. (83). 

We interprete the divergencies of uj at \n\ — > — r — 1 for even r = 0,2, ... as a man- 
ifestation of the double-logarithmic contribution A7 ~ G 2 /uj 3 . Indeed, one can present 
expression (89) in the kinematical region w C 7 C 1 with the same accuracy as a solution 
of the equation 




cj = G - - V(j - 7) + * (1) + G 2 { -U. Al >± + 



2 , - ^(j- 7 ) , c(j-l- 7 ) \ 

7 2 27 ; ' 

where j = 1 + uj + \n\. It is obvious, that if one neglects the term ~ G 2 the anomalous 
dimension 7 contains the double-logarithmic contribution A7 ~ G 2 /uj 3 . Note, that the 
above expression corresponds to a small-7 asymptotics of the BFKL equation written in 
the separable form (cf. (46)) 



1 = c ^ 2*(l)-g( 7 )-g(j-7) , c( 7 -l) + c(j-l- 7 ) | c(2) ^j 
G 2 

+ T (0(7) + 00' -7)) , 

where 0(7) is presented in Eq. (28). On the other hand in the Born approximation at 
j — > — r and 7^0, neglecting the corrections to the relation uj = j + r, one can simplify 
the previous equation as follows 

1 „ q ( 1 _ S 1 (r) + (u - ^) S 2 (r) \ ^ 
\7(^-7) ^ / 

It will be shown below, that for even r the double-logarithmic terms A7 ~ G 2 /uj 3 
obtained from this equation are not cancelled. It is important, that the next-to- leading 
correction ~ G 2 to uj was obtained in ref. [7] by subtracting the first iteration of the Born 
kernel from the contribution of all one-loop diagrams to avoid the double-counting. At 
\n\ -^-r-lwe chose another Born approximation which contains the double-logarithmic 
terms. It means, that now one should subtract from the next-to-leading correction some 
terms including the most singular pole which appears in the first iteration of this Born 
contribution to avoid the double-counting. But the subtracted terms should be small in 
the applicapability region uj <C 7 for the BFKL equation. Thus, the coefficient in front of 
G 2 after this subtraction can be written for even r as follows 

(m-l)- 1 -S 1 (r)-(iu- 1 )S 2 (r) | 2(uj - 7 )' 2 + S 2 (r) + S_ 2 (r) - 2((2) 
uj 7 2 2ury 
1 a(r) b(r) d(r) e(r) 



(uj — 7) 2 7 2 (uj — 7) 7 2 (uj — 7) 2 7 (uj — 7) 7 7 2 ' 

where the residues a(r),b(r),d(r) and e(r) are some unknown functions. The second 
order pole (uj — j)~ 2 and less singular contributions are not essential with our accuracy. 
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Thus, we obtain the following modified BFKL equation at even r and 7 < w 

y'JUJ uj 2 uj J 

l2 fa(r)-S 1 (r) , a(r) + b(r) , 5_ 2 (r) + 3S 2 (r) + C(2) + 2d(r) , e(r)-S 2 (r)\ 



\ 7^0; "yuj z Zuj^y 7^ y 

If one will take into account, that u> and j + r are generally different parameters (see 
(13)), it will be needed to substitute in the above equation uj by j + r and to change the 
unity at its left hand side by an analytic function of j + r 

1 -> 1 + ci(r)(j + r) + c 2 (r)(j + r) 2 + ... , 

where c±(r) = —Ci(r), c 2 (r) = Cf(r) — C* 2 (r) and the parameters Ci(r) and C* 2 (r) are 
defined in Eq. (13). 

After this modification of the equation we calculate the behavior of the anomalous 
dimension at j — > — r for even r in terms of the parameters a(r), 6(r), d(r), e(r) and c 12 (r) 

7 = G (— + tf(r) + L(r) (j + r)^ + G 2 (7—^ + + , (90) 

where 

K(r) = a(r) - S 1 (r) - Cl (r) , T(r) = 6(r) + tf(r) - Cl (r) , 
L(r) = e(r) — S 2 (r) — (a(r) — Si(rfj + c\(r) - c 2 (r) , 

i?(r) = d(r) + ^~ 2(r) + S f r) - C(2) - (a(r) - S^r)) (b(r) + ^(r)) 
- Cl (r)(a(r) + 6(r) - 2S 1 (r)) + 3c 2 (r) - 2c 2 (r) . 

This expression for 7 should be compared with the result (84) obtained in the previous 
subsection from the DGLAP equation. The leading poles Gj (j + r) and G 2 / (j + r) 3 coin- 
cide in both cases. As for the regular terms ~ G and residues of the poles G 2 / (j + r) 2 and 
G 2 /(j + r), the agreement can be achieved with an appropriate choice of the parameters 
a(r) , b(r) , d(r) , e(r) and Ci )2 (r). Thus, to verify the possibility, that the anomalous di- 
mension 7 um (j) can be derived completely from the BFKL equation one should calculate 
these parameters from the Feynman diagrams. As it is seen from the above expressions, 
the parameters Ci ;2 (r) enter together with a(r), b(r), d(r), e(r) and could be omitted. The 
vanishing of c i)2 (r) corresponds to a natural assumption, that to obtain the universal 
anomalous dimension of the twist-2 operators one should initially put \n\ = — r — 1 and 
only after that uj can be pushed to zero. In this case we obtain 7 um (j) (84) if the param- 
eters a(r) , b(r) , d(r) are chosen as follows 

a(r) = b(r) = -S 1 (r + 1), e(r) = ^( r + i)-^2_, 

() _ S 1 (r + 1) g_ 2 (r + l)-3g 2 (r + l) 
1 ' r + 1 2 
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As a result, the equation for the anomalous dimension can be written as follows 

1 = Q ( 1 S 1 (r) + (u - -y) S 2 (r) \ | 
\7(^-7) ^ / 

|C 2/ (r + 1)- 1 fr(r + l) | g 2 (r + 1) - S 2 {r + 1) 

\ 7 2 cj 7 2 (cj — 7) 7(cj — 7) 2 7 2 

| S_ 2 (r) + 3S 2 (r) + ^^ + g_ 2 (r + l)-3g 2 (r + l) + /(r)\ 
2c^7 27(0; — 7) 7 / ' 

where the function /(r) is fixed from the condition, that its solution reproduces correctly 
the regular terms ~ G 2 in 7 um (j) at j — > — r (see (84)). Providing that we shall take into 
account also non-singular terms at 7 — ■> 0, the BFKL equation will contain an important 
information about the singular part of the anomalous dimension in the order G 3 . 

For odd r the situation is simpler because here the double-logarithmic terms at j — > — r 
are absent. Indeed, in this case one can write the BFKL equation in the form (after the 
subsequent substitution uj — > j + r) 

1 + m U + Mr) u 2 = G (— + — ^ — M1±^M ) + G 2 (5W^iW 

\7u; uj 2 uj J \ 7 2 cj 

1 | a(r)+b(r) | -g_ 2 (r) + 3S 2 (r) + 2((2) + 2d(r) | e(r) - ff 2 (r)^ 



uj 3/ y 7W 2 2c^7 7 2 y ' 

where a(r), b(r), d(r), e(r) and c+ 2 (r) are the corresponding parameters for odd r. There- 
fore the behavior of the anomalous dimension at j — > — r for odd r is 

7 = G + *W + l ^ + O) + G 2 (J^y 2 + ^) , (92) 

X(r) = 5(r)-S'i(r)-Ci(r), f(r) = b(r) + 2 K(r) + Ci(r), 

Z(r) = e(r)-5 2 (r)-(a(r)-^(r)) 2 + c 2 (r)-c 2 (r), 

£(r) = d( r ) + e(r) - ^(r) + ^(r) + 2C(2) _ ^ _ + ^ _ 

- ci(r) (a(r) + 6(r) - 2S' 1 (r)) + 25 2 (r) - 5 2 (r) . 

Again we have a qualitative agreement with the result (84) obtained from the DGLAP 
equation, but in the order G 2 the residues of the poles l/(j + r) 2 and l/(j + r) depends 
on parameters S(r), b(r), d(r),e(r) and ci j2 (r). If £i j2 (r) = 0, one can have a complete 
agreement with Eq. (84) provided that 



d(r) = — , b(r) = 2Si(r + 1) , e(r) = S 2 (r + 1) 



2 ' 



r + 1' w 7 ' v ' iV (r + 1) 

d(r) = -5 1 (r)5 1 (r + l) + ^^^ + C(2) + 25 2 (r + l). 

In this case the equation for the anomalous dimension can be written in the form 
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+ (u - i) S 2 (r) \ + 



UJ 



Si(r) (r + 1)- 1 | 2S 1 (r + l) | S 2 (r + 1) - S 2 (r + 1) 
7 2 u; 7 2 (u; — 7) 7CJ 2 7 2 

| 3^ 2 (r)-g_ 2 (r) + 2C(2) | C(2) - S 1 (r)S 1 (r + 1) + j^gg + 2S 2 (r + 1) | /( r )\ 
2c^7 7(0; — 7) 7 / ' 

where /(r) is obtained from the condition that the regular terms A7 ~ G 2 at j — > — r 
coincide with the corresponding contributions to A^y um . Taking into account also the non- 
singular contributions in the BFKL equation at 7 — > one can obtain some information 
about singularities of 7 in the order G 3 . 

Thus, the results for the universal anomalous dimension r ) uni {j) obtained from the 
BFKL equation with the use of our hypothesis, that r y um (j) coincides with r ) BFKL {j) af- 
ter its analytic continuation to the singular points \n\ = — r — 1, agree in main features 
with the expressions derived directly from the DGLAP equation, but for a full agreement 
one should verify by calculating the corresponding Feynman diagrams, that the parame- 
ters a(r), b(r), d(r), ci j2 (r) and 5(r), b(r), d(r), ci j2 (r) reproduce correctly the non-leading 
singularities G 2 / (j + r) 2 and G 2 / (j + r). 



7 Conclusion 

Above we reviewed the LLA results for the anomalous dimensions of twist-2 operators in 
the N=4 supersymmetric gauge theory and constructed the operators with a multiplica- 
tive renormalization (cf. [28]). These anomalous dimensions can be obtained from the 
universal anomalous dimension r y um {j) by a shift of its argument j — > j + k because they 
belong to the same supermultiplet. We calculated the NLO correction to 7 um (j) by using 
a number of plausible arguments. The results have a compact form (see eqs. (70) and 
(68)) in terms of polylogarithms and associated special functions. They are verified by 
direct calculations [42] . Note that recently the LLA anomalous dimensions in this theory 
were constructed in ref. [54] in the limit j — > 00 from the superstring model with the use 
of the AdS/CFT correspondence [55, 56, 57]. 

We investigated properties of the next-to-leading corrections to the kernel of the BFKL 
equation in the N=4 supersymmetric theory. The absence of the coupling constant renor- 
malization in this model leads presumable to the Mobius invariance of the BFKL equation 
in higher orders of the perturbation theory. But the holomorphic separability of the BFKL 
kernel is violated in the NLO approximation (see (BIO)). Instead we have the hermitial 
separability of the corresponding Bethe-Salpeter kernel (see (46)). The cancellation of 
non-analytic contributions proportional to 5® and <5 2 in = 4 SUSY is remarkable (such 
terms contribute to uj in QCD and in N=l,2 supersymmetric models [9, 53]). Moreover, 
the hermitian separability of the BFKL kernel is also valid only in N=4 SUSY. 

These properties could be a possible manifestation of the integrability of the reggeon 
dynamics in the Maldacena model [55] (see also [56, 57]) which is reduced to the classical 
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supergravity in the limit N c — > oo. Indeed, in this N=4 supersymmetric model the 
eigenvalues of the LLA pair kernels in the evolution equation for the matrix elements 
of the quasi-partonic operators are proportional to ^(j — 1) — [28], which means, 
that the corresponding Hamiltonian at large N c coincides with the local Hamiltonian 
for an integrable Heisenberg spin model [19]. The residues of these eigenvalues at the 
points j = —r were obtained in Ref. [9] from the BFKL equation in LLA by an analytic 
continuation of the anomalous dimensions to negative integer values of the conformal 
spin |n|. It was shown above, that the analogous correspondence between the BFKL 
and DGLAP equations takes place qualitatively in the NLO approximation (cf. Eq. 
(84) and Eqs. (90, 92)). To verify quantitatively in this approximation the hypothesis, 
that in N=4 SUSY the DGLAP equation can be obtained from the BFKL equation, one 
should calculate the next-to-leading corrections to the BFKL equation in the region, where 
|n| + r + 1 — > and 7 ~ u. We hope to return to this problem in our future publications. 
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8 Appendix A 

Here we consider the most complicated contribution $(|n|,7) + $(|n|, 1 — 7) to the NLO 
correction 8(n, 7) in the BFKL equation (21) and derive the property similar to the 
holomorphic separability (cf. [14]). To begin with, let us split the function $(n, 7) in two 
terms 



$(|n|, 7 ) = $ 1 (|n|,7) + $ 2 (|n|,7), 

where 

00 (P'(k + \n\ + 1) - (-l) k m\k + \n\ + 1)) 

*><H,7)= E ftm L 

°° (-l) fc (^(|n| + k + 1) - 

and 
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~ + 1) + + 1)) 



"go V + M) 2 

It is convenient also to write the functions <3>i(|n|,7) and <£> 2 (|n|,7) as 



$i(|n|,7) = $ (1) (M,M) - $ (2) (|n|,M) , 
$ 2 (N,7) = $ (1) (0,M) + $( 2 )(0,M), 



where 



*V(\n\,s) = £ 



+ \n\ + 1) 



fc=0 

OO / 1 \fc 

$( 2 )(|n|,s) = £ i-^L ( tt'(Jfe + |n| + 



fc=0 



_ ^(H + fc + i)-^(i) \ 



(Al) 

(A2) 
(A3) 



If we expand fi'(k + \n\ + 1) in the series 

OO 

£(-l) m /(' + A:+M + l) 2 
z=o 

and insert it in the sum over k in the r.h.s of Eq.(A2), we obtain 

$ (1) (M,s) = $ (2) (M, \n\ - s + 1) + -s + l)[*(l) - *(s) 
In particular, 

$ {1) (|n|,M) = $ (2) (|n|,l - M) + - M) [tf (1) - (M)], 

$ (1) (0,M) = $ (2) (0, 1-M) + - M) - *(M)" . 

The function <3> 2 (|n|,7) depends only on M and therefore the corresponding contri- 
bution to uj has the property of the hermitian separability. Further, for $i(|n|,7) and 
( &2(|^|,7) we derive the following relation 

$i(|n|,7) + $i(|n|,l-7) 

= M) - $ (2) (|n|, M) + 1-M) - $ (2) (M, 1 - M) 

= /?'(M) [*(1) - *(1 - M)] + - M) [*(1) - *(M)] , (A4) 



$ 2 (M) = $ (1) (0,M) + $ (1) (0, l-M)-ff{M)\V{l)-V(l-M)\ 
= $ (2) (0, M) + $ (2) (0, 1 — M) + - M) [tf (1) - *(M)] 



(A5) 
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Therefore one can obtain 

$(|n|, 7) + $(M, 1 - 7) = X(n, 7) (P'(M) + f3'(l - M)) 

+$ 2 (M) - (3'(M) [tf (1) - *(M)] + $ 2 (1 — M) — - M) [tf (1) - - M) 
where x(. n il) is given by Eq.(20). 



9 Appendix B 



In this Appendix we demonstrate a violation of the generalized holomorhical separability 
in the NLO correction 5(71,7), i- e - we show an impossibility to present Eqs.(25)-(28) in 
the separable form symmetric to the substitution m (i.e., respectively, M <-> M for 
both positive and negative \n\). 

The expressions appearing in S(n, 7) (21) can be written as follows 



*"(M) + - M) 



1 



(Bl) 



*"(M) + — M) + *"(M) + - M) 

$(|n|, 7 ) + $(|n|, 1 - 7) = $ (2) (0, M) + $ (2) (0, 1 — M) + $ (2) (0, M) + $ (2) (0, 1 - M) 
+ (/3'(l — M) + - M)) [tf (1) - *(M)] + (p'(M) + (3'{M)) [tf (1) - (1 - M)J .(B2) 

Using the property 



sin 2 (Mvr) ' 

we can present the terms containing the functions (3' in the following form 
(3'(1 — M) + (3'(1 — M) = 



cos(M7r) 



— 7T 



f3'(l — M) + p'(l — M) + P'(M) + f3\M) + 7T 2 - 

P'(M)+p'(M) = 
(3'{M) + (3'(M) + (3'{l — M) + (3'{1 — M) — tt 2 C ° S ^\ + tt 2 



cos((l - M)tt) 



sin 2 (M7r) " sin 2 ((l-M)7r) 



cos((l - M)?r) 



sin 2 (M7r) S in 2 ((l - M)vr) 



Then, the second line in Eq.(B2) can be replaced by 
1 



2 
+ 



(5'{M) + (5'{M) + — M) + /3'(1 - M) 

1 /7T 2 COs(M7r) 7T 2 C0S((1 - M)tt; 



2 V sin 2 (M7r) s i n 2 ((l - M)tt) 
where x( n ij) from (19) has the symmetric representation: 

1 



— M) — *(M) 



X (n, 7 ) = 2tt(l) 



*(M) + (1 — M) + *(M) + - M) 



(B3) 



(B4) 
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Thus, the NLO term S(n, 7) is presented as follows 



6(n,i) = «$i(n,7) + £ 2 (ra,7), 



where 



6 1 (n n ) = <f>(M)+<f>(M)-^(p(M)+p(M)) ("0 = 4a X (n, 7)), 

/ 7T 2 COs(M7r) 7T 2 C0S((1 - M)^ 



^2(^,7) 



sin 2 (M7r) sin 2 ((l - M)ir) 



*(M)-*(1-M) 



and 



2p(M) 
20(M) 



= P'(M) + — M) + ((2) , 



(B5) 
(B6) 

(B7) 



= 6C(3) + *"(M) + *"(1-M)-2$ (2) (0,M)-2$ (2) (0,1-M). (B8) 



All above terms are symmetric to the substitution M «-> (1 — M) in an agreement 
with previous results. Moreover, the term Si(n, 7) is symmetric to M <-> M (and, thus, 
to n <-> — n) due to the representations (B7) and (B8). Therefore for this contribution 
the property of generalized holomorphic separability is valid. 

The term 5 2 (n, 7) can be presented as 



/ 7T 2 COs(M7r) 7T 2 COs(M7t) 



sin 2 (M7r) sin 2 (M7r) 



*(M) - *(M) 



7T COs(M7r) 



sin(M7r) 



where 



7TCOs(M7r) 7TCOs(M7r) 



sin(M7r) sin(MTr) 
It can be splitted in two parts: S 2 (n, 7) = S^\n, 7) + S^\n, 7) where the contribution 



#>(n, 7 ) 



7T 2 COs(M7r) 7T 2 COs(M7r) \ 7T COs(M7t) 



sin 2 (Mvr) ' s in 2 (Mvr) / sin(MTr) 

7T 3 COS 2 (M7r) 7T 3 COS 2 (Mn) 



siu^Mvr) 



sin 3 (M7r) 



is symmetric under the substitution M <-> M (and, thus, to n <-> — n). So, the term 
#2^ (h, 7) has also the property of the generalized holomorphic separability. 

Thus, the sum of the contributions Si (n, 7) and S^ (n, 7) has the generalized separable 
form symmetric to the substitution m <-> m: 



Si (n, 7) + S ( 2 ] (n, 7) = (f)(m) + </>(m) - ^ ( p{m) + p(m) 



— 71" 



'cos 2 (m7r) cos 2 (m7r) 



+ 



sin (mn) sin (run) 



(B9) 
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The last term 



\ sm (Mir) 



+ 



7T 2 C0s(M7r) 

sin 2 (M7r) 



*(M) - *(M) 



is anti-symmetric under the transformation M <-> M, i.e. 
when n > 0, and thus, M = m, M — m it equals to 



(7r 2 cos(m7r) 7r 2 cos(m7r) 
sin 2 (m7r) sin 2 (m7r) 

and when n < 0, and thus, M = m, M = 
( 7r 2 cos(m7r) 7T 2 cosfm7r" 



sm (mir) 



+ 



sin (m7r) 



m it is 

\&(ra) - *(m) 



Therefore the term 5^ (n, 7) in the next-to-leading correction 



<5(n,7) 



fm) + 



m 



2a 



p(m) + p(m) 



7T 



cos 2 (m7r) cos 2 (m7r) 
sin 3 (m7r) sin 3 (m7r) 



+ *? ) (n,7) 



(BIO) 



violates the holomorphic separability. Note, that the anomalous term 5^\n,^f) can be 
written as Eq. (31), i.e. it is zero for odd n, where 8(n, 7) coincides with its analytic 
continuation to corresponding negative \n\. 



10 Appendix C 

Here we construct the anomalous dimension matrices 7^0) ( a > b = g,q,tp) and 7^(7) 
(a, b — g,q) in terms of their eigenvalues in LLA 

7 i 0) (j) = -45x0' =F 2), 7^0') = -^(j), 7i%') = -4Si(j =F 1) • 
To illustrate the procedure we consider the polarized and unpolarized cases separately. 

Polarized case. From Eqs. (58), (59) and (61) we have 



TSHj) + ^(J) = if^j) + 27^0') = 7fO') = -4S!(j - 1) (CI) 



T^OV^y^a) = T^OO - 7i?(i) = 7i 0) (j)--4^(j + l) (C2) 

It is obvious, that Eqs.(55)-(59) and (C1)-(C2) correspond to a diagonalization of the 
anomalous dimension matrix. So we can rewrite Eqs. (C1)-(C2) as 

' tSPW) t?W) ' €»0) l l ' ^ v • (C3) 
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where 



V = 



v g -2 

\ V 9 



3-1 



j+2 J 



2j + l\-\v- 1 v- 1 



(C4) 



with arbitrary values of v g and v q . 

Thus, Eq.(C3) leads to the following representation of the anomalous dimension matrix 
in the polarized case 



~(o),., 
7 0) 



T$2U) t|?0) 



T^O) 



V 









7< 0) 0) 



V 



7f(i) 



2(j + 2) 4(j - 1) 



+ 7<°>0) 



2(j" - 1) -40' -1) 



2(2j + l) L' + W/ V ^' + 2 20-1); '" wy V-0' + 2) 20 + 2) 
which does not depend on v g and v q . 

Unpolarized case. From Eqs. (58), (59) and (60) we have 



(C5) 



7ff 0) + 7g ; 0) + 7$ 0) = 7g' 0) + 7ff 0) + 7S? 0) 
= 7$0)+7«?0)+7<!?0) = 7f0) = 



45i0-2), 



(C6) 



1 



J + 1 



7$0) 



! 99 ^ J > 2(j — 1) qg 3(j — 1) l<P9 



= 7i?0) - 20 - l)7i?0) + ^^7^0) 

7$0) = 7S 0) 0)^-4^iO), 



7^ 0) 7 + 1 ^¥>(^') + 



3 
3 



20 + 1) 



(C7) 



7 



SO) 



J+ %;^+n^7so) 



= 7j?0) 



- 1) q9 

3-1 



J + 2 



7 



S?0) 



j'0 - 1) 



7^0) + 7-^ ^ 0)- 



-7S0) 



7 i° ) 0) = -45 1 + 2). (C8) 



+ 1)0 + 2)'"^' j + 
The formulae (C6)-(C8) are equivalent to the matrix equation 

I \ 





7fi>0) 7^0) 7^0) 
V^ 1 I 7|0) 7|0) 7<°>0) | V 

7g<^ ) 7<jr< ( j ) 7(^(^5 ) 



V 







7$ 0) 0) 
7 ( ° ) 0) J 



(V- 1 V = l) ,(C9) 



where 



V 



-20 " 1>. 



\ V 9 |0 + 



(CIO) 
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V 



-1 



2(4j'2-1)(2j + 3) 



/ (2j + 3)v-i 4^(2j + 3W- 1 



6 j+l V q 



i+2 

Q_M±1 



0+1)0+2) 



V (2j - l)^ 1 4(2j - l)v, 



-1 



j(2j+l) ..-1 
0+1)0+2) U 9 

3(2j - 1)< 



with arbitrary values of i> g , t> g and i^. 

Analogously to the polarized case one can write the expression for the anomalous 
dimension matrix in the form 



7!! } (j) 

0)1 



7<?(j) 



\ 



I if ij) 
if (3) 

7 W (J) J 





-,(«)/ 



2(4j'2-1)(2j + 3) 



+ (2j + l)lf(j) 



(2j + 3) 7 f(j) 



/(j + 1)(j + 2) 4(j 2 
(j + l)(j + 2) 4(j 2 
V (j + 1)(j + 2) 4(j 2 



/ 6(i-l)0' + 2) -12(j-l) -6j(j 
-3(j + 2) 6 3j 

-2(j + l)0' + 2) 4(j + l) 2j(j + l 



1) 
1) 
1) 

1) \ 



3J'(J " 1) 
3j(j - 1) 
3j(j - 1) 



/ 



+ (2j-l) 7 i 0) (j)- 
which does not depend on v g , v q and v v . 



3(3 - 1) !./(./ - 1) 3j(j - 1) 

-3(3 + 2) -43(3 + 2) 33(3 + 2) 

(j + l)(j + 2) 4(j + l)(j + 2) 3(j + 1)0' + 2) 



(Cll) 



11 Appendix D 

Here we obtain the expression (16) for the t-channel partial wave f w (n) of the process 
e + e~ — > fi + fi~ in the double-logarithmic approximation. Note that our results coincide 
with those obtained in the original paper [27]. 
The Bethe-Salpeter equation (15) 

r s tik'? /•mm(s,sk'?/k 2 . ) A J / rv \ 

can be transformed after going in the tu-representation 

Ja-ioo 2niu \kj_J \ m 2 J 

to the form 

/«(«) = l + ^(jy^'e-^ K '-^UK') + j^dK'UK')\ . 
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We introduced here the ultraviolet cut-off A in lnk±. It will be pushed to infinity in 
the end of calculations. Note, that in QCD such cut-off could be interpreted as an effect 
of the interaction vanishing at large k 2 due to the asymptotic freedom. The differentiation 
of the integral equation gives 



d_ 



/*(«) = A ^ dn' e~^' U{k') . 

J K 



The second differentiation leads to the differential equation 

We can find its solution in terms of two arbitrary constants 
where 



(Dl) 



/From the above expression for e u}K -j-f ul (K) we obtain, that 



and from the integral equation 



1 = 



A / a 



Therefore 



n+ -- p A ^ 

A 



lo> 101 



1 -1 



> a u 



A 



7^ 
loi 



-i -1 



In particular, for A — > oo we obtain 

where the t-channel partial wave on the mass shell 



satisfies the infrared evolution equation 

A 



/« = l + ^/«j- (D3) 
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In this case the partial wave can be written also as follows 



f%+i°° cfy / k 2 



) 



7 



(u - 27)7^ X' 2 




(D4) 



Note, that providing that A = oo from the beginning the condition a+ = follows 
from the requirement, that the t-channel partial wave f w {n) does not grow at cu — > oo (in 
an opposite case we can not reduce the initial Bethe-Salpeter equation for A(s, k\) to the 
equation for because the lower limit of integration s' = k' 2 is not zero). 

Thus, we derived the correct expression (16) for / w (7). The '— ' component of the 
value (Dl) coincides with one (17) in the section 2. 
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